Velocity Profiles in Slowly Sheared Bubble Rafts 
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Measurements of average velocity profiles in a bubble raft subjected to slow, steady-shear demon- 
strate the coexistence between a flowing state and a jammed state similar to that observed for 
three-dimensional foams and emulsions [Coussot et al„ Phys. Rev. Lett. 88, 218301 (2002)]. For 
sufficiently slow shear, the flow is generated by nonlinear topological rearrangements. We report on 
the connection between this short-time motion of the bubbles and the long-time averages. We find 
that velocity profiles for individual rearrangement events fluctuate, but a smooth, average velocity 
is reached after averaging over only a relatively few events. 

PACS numbers: 83.80.Iz,83.60.La,64.70.Dv 
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A proposed jamming phase diagram which treats 
applied stress in a manner analogous to temperature and 
density, offers an interesting framework for the study of 
a wide range of systems subjected to shear 0. The jam- 
ming phase diagram proposes the existence of a new state 
of matter, the "jammed state", and suggests that similar 
transitions to this state would occur as a function of ap- 
plied stress, temperature, or density. For example, many 
complex fluids, such as emulsions, foams, and granular 
materials, exhibit a yield stress below which the mate- 
rial does not flow, or jams. As one approaches the yield 
stress from above, the viscosity diverges. This is analo- 
gous to the behavior of the viscosity in the glass transi- 
tion. Experiments have confirmed the applicability of the 
jamming phase diagram in certain colloidal systems Q. 
These are examples of continuous jamming transitions. 

Recent work in a range of soft materials, including 
colloids, granular matter, foams, and emulsions, suggest 
that the jamming transition can be discontinuous • 
The materials exhibited a yield stress. Yet, under con- 
ditions of steady shear flow, there existed a critical ra- 
dius at which the shear-rate was discontinuous and a 
transition from a flowing to a jammed state occurred 
0. This is in contrast to the results of most stress 
versus rate of strain measurements for materials with a 
yield stress. Such experiments suggest the stress is well- 
modelled as a continuous function of shear-rate, such as 
with a Herschel-Bulkley model |6j. Additional evidence 
for the discontinuous transition is given by the observa- 
tion of a bifurcation in the material's viscosity during 
constant stress experiments 0. 

In addition to the connection with the jamming tran- 
sition, the coexistence of a flowing and jammed state is 
one example of another feature of some jammed systems: 
shear localization 0, B IS Ell EH E3 • In general, shear 
localization (or banding) refers to a flow state which is 
spatially separated into a regions of different shear rates. 
For wormlike micelles, shear banding has been observed 
in which there is a discontinuous transition from a high 



shear rate region to a low shear rate region |l2( ■ Experi- 
ments in granular systems 0, H, 0, E| and confined foam 
[To| report a type of shear localization in which the veloc- 
ity profile is exponential. In these cases, the shear rate is 
continuous, though the system is effectively divided into 
high shear rate and zero shear rate regions. Simulations 
of Lennard- Jones particles and quasi-static foam 01 
confirm that shear localization can occur independent of 
the existence of a yield stress. In the case of quasi-static 
foam, the simulations suggest that a key element is the 
localization of slip events ^ij. Understanding the dif- 
ferent sources of shear-localization will contribute to a 
better understanding of the jamming phase diagram. 

In this Letter, we confirm the coexistence between a 
flowing and a jammed state for slow, steady shear of a 
bubble raft. By considering the behavior on relatively 
short time scales, we show the connection between fluc- 
tuations that occur in the stress and the nonlinear flow 
events that comprise the unjammed state. This provides 
insight into the connection between the individual bubble 
motions and the observed average flow properties. The 
average stress versus rate of strain curves have been mea- 
sured separately These curves are consistent with a 
Herschel-Bulkley model for the system. However, in light 
of the velocity measurements reported here, it is impor- 
tant to reconsider the interpretation of those results. 

A bubble raft consists of a single layer of bubbles float- 
ing on a fluid surface ^(|. Our bubble raft is described 
in detail in Ref. 0|. A random distribution of bubble 
sizes was used, with an average radius of 1 mm. The 
Couette viscometer is described in detail in Ref. . To 
shear the foam, an outer Teflon barrier is rotated at a 
constant angular velocity. We measure the azimuthal ve- 
locity, vg(r). The shear rate is given by j(r) — r-^^^- 
Therefore, we considered the normalized angular ve- 
locity, v(r) = vg(r)/(ilr), where the jammed state is 
v(r) = 1, or 7 = 0. The stress, cr(rj), on the inner 
cylinder was monitored by measuring the torque, T, on 
the inner cylinder: <r(rj) = T/(2irrf). The inner cylin- 
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der was supported by a torsion wire (torsion constant 
k = 570 dyne cm), and T was determined from the an- 
gular position of the inner cylinder. Therefore, the in- 
ner cylinder had an instantaneous angular speed, even 
though its average angular speed is zero. At both bound- 
aries, the first layer of bubbles was never observed to slip 
relative to the boundary. This feature combined with the 
finite size of bubbles set the effective inner (r,; = 4.3 cm) 
and outer (R = 7.2 cm) radii. We report results for two 
rotation rates of the outer cylinder, £1 = 8 x 10 -4 rad/s 
and f2 = 5 x 10 -3 rad/s, which corresponds to shear rates 
(at r = 4.3 cm): 7 = 4 x 1(T 3 s" 1 and 7 = 3 x 1(T 2 s" 1 . 

The < <r > versus 7 behavior was measured separately 
and is reported elsewhere It is consistent with a 

Herschel-Bulkley form for the viscosity: < cr(ri) >= 
Hoi n + t , with n = 1/3, and r Q = 0.8 ± 0.1 dyne/cm. 
For 7 < 0.1 s -1 , the < a > was essentially independent 
of shear rate, though it fluctuated between 0.5 dyne/cm 
and 2 dyne/cm |l5j. Both shear rates reported on here 
are within this quasi-static limit. The gas area- fraction 
was 0.95. As the bubbles are actually three-dimensional, 
we used an operational definition of gas area-fraction: 
the ratio of area inside bubbles to the total area of the 
bubble raft in a digitized image. 

The fluid substrate (subphase) is driven at the same 
time as the bubbles [lj|. Tests were made with bubble 
rafts that did not touch the outer barrier. Under rota- 
tion of the outer barrier, no motion of the bubbles was 
detected, ruling out driving by the subphase. In addi- 
tion, because the subphase is sheared with the bubbles, 
the bubbles and subphase have similar velocities during 
most of the motion. This minimizes any viscous dissi- 
pation between bubbles and water. Also, the effective 
internal viscosity of the bubble raft is 300 to 2500 g/cm s 
for the range of 7 studied. This is a factor of 10 4 to 
10 5 greater than the viscosity of the subphase. This ra- 
tio ensures that dissipation between bubbles dominates 
dissipation from the subphase-bubble interactions. 

Bubble velocities were measured from taped video im- 
ages that were digitized in the computer. For ft — 

8 x 10~ 4 rad/s, the time between digitized images was 
3.2 s/image, and for — 5 x 10~ 3 rad/s, it was 2.0 
s/image. An image processing routine based on standard 
National Instruments Labwindows™/CVI functions de- 
tected and tracked individual bubbles. The velocity was 
calculated using the total displacement of 10 digitized im- 
ages. This allowed tracking of the rapidly moving bub- 
bles from image to image (short time between individ- 
ual images) and sufficient total time to measure bubbles 
as slow as 6 x 10~ 4 cm/s for = 8 x 10~ 4 rad/s and 

9 x 10~ 4 cm/s for = 5 x 10~ 3 rad/s. Finally, the bub- 
bles within each of 24 equally spaced radial bands in the 
range 4.3 < r < 7.2 cm were averaged to compute v(r). 

Figure H is a plot of a(ri) versus time for the velocity 
data reported on here. The initial elastic regime consists 
of a linear increase of a(ri) with time. The subsequent 
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FIG. 1: Stress versus time for f2 = 5 x 10 -4 rad/s [gray line 
labelled (1)] and Q. = 8 x rad/s [black line labelled (2)]. 
The solid horizontal line is t determined from a fit using 
data from Ref. 1.5J to a Herschel-Bulkley fluid model. The 
dashed horizontal line is a(r^) when cr(6.7 cm) = t (black 
line) or er(6.3 cm) = r (gray line). The dotted line is a(ri) 
when a(R) = r . The vertical lines indicate the stress drops 
presented in Fig.|3]for Q = 8 x 10 -4 rad/s. 



"flowing" regime is dominated by irregular variations in 
the stress characteristic of the slow shear-rate "flow" of 
many jammed systems |2J. For comparison with the flow 
data presented later, the solid line in Fig. ^ represents 
t (the "yield stress") based on the fit of < u(ri) > to 
a Herschel-Bulkley model ^ij. (Given the existence of a 
discontinuous transition, such a definition of a yield stress 
may not be meaningful.) For each data set, the dashed 
line represents the value of cr(ri) such that o~(r c ) — t d , 
assuming a(r) = T/(2irr 2 ) [|J. (r c is the radius at which 
the system jams, see Fig. The dotted line is cr(ri) 
such that a(R) = t d . 

Figure |21 shows v (r) versus r for fl — 5 x 10~ 3 rad/s for 
a number of different averages. In each case, the average 
is computed using bubbles from roughly 1/3 of the sys- 
tem. The solid circles represent an average over 1000 s, 
starting 210 s after the initiation of shear. This corre- 
sponds to approximately 2800 individual velocities per 
channel. Defining an event as a consecutive period of 
stress increase and decrease, the other curves are aver- 
aged over a single event (□), four events (o), 10 events 
(A), and 20 events (V), respectively. The 10 event aver- 
age is in reasonable agreement with 1000 s average, and 
the 20 event average (400 s) is indistinguishable from the 
1000 s average. Figure [5] also shows v(r) versus r for 
il = 8 x 10 -4 rad/s. The average covered a total time of 
1020 s, starting 650 s after the initiation of shear. This 
corresponds to approximately 2000 individual velocities 
per channel and on the order of 50 events. 

To find 7 and r c , the average velocity is fit to v(r) = 
A + B/r 2 / n (solid curves) over the range < v(r) < 
0.98. This is the velocity profile for a power-law fluid 
(cr(r) oc r") in a Couette geometry 0. A number of 



3 



T 1 1 1 1 1 1 — 1 




radial position (cm) 

FIG. 2: v(r) = vg(r)/(rSl) versus r for SI = 5 x 10~ 3 rad/s 
averaged over time from t — 250 s to the end of the run (•) 
and for SI = 8 x 10 -4 rad/s averaged over time from t = 650 s 
to the end of the run (■). The solid lines are fits to a power- 
law model for viscosity. The dashed line is the v = 1 line. 
Also shown are v(r) for SI = 5 x 10~ 3 rad/s averaged over a 
single event (□), four events (o), 10 events (a), and 20 events 
(V). The insert illustrates the discontinuity in the shear rate 
for SI = 5 x 10~ 3 rad/s (•), SI = 8 x 10~ 4 rad/s (a), and 
SI = 1 x 10~ 4 rad/s (□). 



aspects of this fit differ from the expected solution for a 
Herschel-Bulkley fluid. First, there is a discontinuity in 
7. This is highlighted by the insert in Fig. [3 Here data 
is shown for three rotation rates, SI = 5 x 10~ 3 rad/s, 
!! = 8x 1CT 4 rad/s, and ft = 1 x 1CT 4 rad/s. (For 
ft = 1 x 10~ 4 rad/s, there are only approximately 10 
events. Based on the results for SI — 5 x 1CP 3 rad/s, this 
is sufficient to confirm the discontinuity.) Quantitatively, 
the crossing of the fit to a power-law velocity profile and 
v(r) = 1 defines the critical radius, r c , and critical shear 
rate, j(r c ). For SI = 8 X 10~ 4 rad/s, r c — 6.7 cm, and 
7(r c ) = 6 x 1(T 4 s _1 . For ft = 5 x 1(T 3 rad/s, r c 
()..'! cm, and 7 fry) = 4 x 10~ 3 s . Unlike in some foam 
experiments |ljj, the critical shear rate differs for the 
two speeds. Second, the fits give n = 0.45 ± 0.05 for ft = 
8 x 10" 4 rad/s and n = 0.33±0.02 for ft = 5 x 10" 3 rad/s. 
For the faster rotation rate, the exponent is in agreement 
with the exponent in the Herschel-Bulkley fit to the stress 
W| and velocity profiles measured at higher shear rates 
17]. However, the measured exponent is different for the 
two rotation rates, as seen in other systems Q. 

The trend for average r c is opposite the expectation 
for a Herschel-Bulkley fluid, where a continuous decrease 
in < a > results in a continuous decrease in r c . It is 
consistent with the measured behavior of < a(ri) > in 
the quasi-static regime. For ft = 8 x 10~ 4 rad/s, the 
system spends more time at a higher value of stress (see 
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FIG. 3: Average velocity during the stress drops indicated 
in Fig. for SI = 8 x 10 -4 rad/s. The inserts are images of a 
subset of the bubbles, with gray circles corresponding to v(r) 
in the same sense as SI, black circles corresponding to v(r) 
opposite SI, and white circles \v(r)\ < 6 x 10 -4 cm/s. 



Fig.^J, and r c is larger. This is also different then trends 
reported in Ref. |j] and [l9(- for higher shear-rates in 
the "continuum" regime. This difference is not surpris- 
ing given that the behavior in the continuum regime is 
"smoother" than in the quasi-static, or "discrete" , regime 
in which we worked. 

Figure [3] highlights the importance of considering 
short-time velocity profiles, in addition to the average 
quantities. In Fig. |3| plots of velocity profiles averaged 
over a single stress drop and corresponding snapshots of 
the bubble motions are shown. The stress drops corre- 
sponding to Figs.OJa)-(d) are indicated by vertical lines 
in Fig. ^ labelled with corresponding letters. Here, the 
shear-rate discontinuity is more apparent. The individ- 
ual velocity profiles are highly nonlinear and not consis- 
tent with a simple continuum model for viscosity. As 
expected r c (the transition point to elastic flow) fluctu- 
ates. However, these fluctuations are not consistent with 
the continuum expectation of a oc 1/r 2 that predicts 
r c = [T/(2ttt )] 1 / 2 , with T = 27rr 2 cr(r l ) 6]. For such 
a model, stress drop (C) would have the largest value 
of r c given its value of er(rj). However, r c is greater for 
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both (D) and (B). This behavior is indicative of stress 
chains, or other nonuniform stress distributions, existing 
in the foam, similar to those observed for granular disks 
in two-dimensions 0]. 

Snapshots of the selected bubble motions are presented 
as inserts in Fig. The images are color coded so that 
bubbles moving opposite the outer cylinder are black, 
gray bubbles are moving in the direction of rotation, and 
white bubbles have \v(r)\ < 6 x 1CP 4 cm/s. The row of 
bubbles at each boundary is not shown. The snapshots 
confirm that r c is connected to bubble rearrangements. 
Both (B) and (D) have the largest r c and radial posi- 
tion at which negative velocities are observed (see circled 
region in Fig. Ot). 

In summary, we have observed the coexistence of a 
jammed and a flowing state in a bubble raft in the quasi- 
static shear limit, as has been observed for other soft 
matter systems at higher shear-rates Q. For the aver- 
age velocity profiles, the transition appears to be discon- 
tinuous and occurs at a critical radius set by the yield 
stress. These are two ways that the average profiles for 
the bubble raft differ from observed exponential shear- 
localization in confined, two-dimensional foam and 
granular systems 0, d- F° r these systems, the 
critical radius is significantly smaller and the shear-rate is 
continuous. At this point, the reasons for the differences 
are not clear, though one potential candidate is the role 
of viscous dissipation. For the confined foam, simulations 
suggest that the exponential velocity profile observed is 
due to localization of the nonlinear rearrangements L4|. 
These simulations do not include viscous dissipation 14| . 
Because this localization is not observed in the bubble 
raft (see Fig. |31 especially event C and D), there may 
exist differences in the role of viscous dissipation in the 
bubble raft and the confined foam. Another potential 
difference is the yield strain, which is quite large in the 
bubble raft system ^?|]- This could impact the distribu- 
tion of rearrangement events. The results of Fig. point 
to the need to understand the velocity profiles during 
individual events in order to resolve these outstanding 
issues regarding the average behavior. 

The results presented in Fig. EH raise two important 
questions. First, given the highly nonlinear and fluctu- 
ating character of the velocities during individual events 
(see Fig. |3J), why does the average velocity converge to 
a smooth curve after averaging over only 20 such events 
(and even fewer)? Second, what sets the critical radius 
for the individual events (Fig. and how is the distri- 
bution of r c for these events related to the value of r c 
found for the long-time averages (Fig. The deter- 
mination of the critical radius is particularly interesting 
given the results presented in Fig. [3] and the work in 
other systems. Unlike the work in the continuum limit 
, r c does not appear to be set by a critical shear rate ei- 
ther for the average profiles (see Fig.[2J or the short-time 
profiles. In fact, the concept of shear rate is not well- 



defined for the motion during individual events, and yet 
each event has a well-defined value of r c . For individual 
events, more work is needed to understand the connec- 
tion between stress and r c . However, r c is clearly not set 
by the stress on the inner boundary, suggesting the need 
to understand the stress distribution within the bubble 
raft. Finally, understanding the connection between the 
individual and average quantities is necessary to resolve 
the apparent discrepancy between < a > versus 7 mea- 
surements that suggest a continuous transition and the 
discontinuous velocity profiles presented here. This is 
necessary for deeper insight into the implications of the 
jamming phase diagram for slow, steady-shear. 
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